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We analyze a microscopic RNA model, which includes two widely used models as limiting cases, 
namely it contains terms for bond as well as for stacking energies. We numerically investigate 
possible changes in the qualitative and quantitative behaviour while going from one model to the 
other; in particular we test, whether a transition occurs, when continuously moving from one model 
to the other. For this we calculate various thermodynamic quantities, both at zero temperature 
as well as at finite temperatures. All calculations can be done efficiently in polynomial time by a 
dynamic programming algorithm. We do not find a sign for transition between the models, but 
the critical exponent of the correlation length, describing the phase transition in all models to an 
ordered low-temperature phase, seems to depend continuously on the model. Finally, we apply the 
e-coupling method, to study low excitations. The exponent 6 describing the energy-scaling of the 
excitations seems to depend not much on the energy model. 



PACS numbers: 64.60.Fr, 87.15.Aa 
I. INTRODUCTION 

RNA plays an important rule in the biochemistry of 
all living systems (2]. Similar to the DNA, it is a lin- 
ear chain- molecule build from four types of bases, i.e., 
adenine (A), cytosine (C), guanine (G), and uracil (U). 
It does not only transmit pure genetic information, but, 
e.g., works as a catalyst. While for the former the pri- 
mary structure, i.e., the sequence of the bases, is relevant, 
for the later the kind of higher order structures, i.e., sec- 
ondary and tertiary structures, are relevant. 

Like in the double helix of the DNA, in RNA comple- 
mentary bases can build hydrogen bonds between each 
other. Compared to DNA, where the bonds are built 
between two different strands, RNA builds bonds be- 
tween bases of the same RNA strand. The information, 
which bases of the strand are paired, gives the secondary 
structure, and the spatial structure is called the tertiary 
structure. The tertiary structure is stabilized by a much 
weaker interaction than the secondary structure. This 
leads to a separation of energy scales between secondary 
and tertiary structure, and gives the justification to ne- 
glect the later 0]. Therefore we deal here with the sec- 
ondary structure only. 

One crucial point in calculating the secondary struc- 
ture is the used energy model: On the one hand, if one 
aims to get minimum structures close to the experimen- 
tally observed one, one uses energy models that take into 
account structure elements 0,|MEII3> ^-S-' ^i^i^' loops. 
On the other hand, if one is interested in the qualitative 
behaviour, one uses models as simple as possible that 
keep the general behaviour, e.g., only one kind of base 
or using energies depending only on the number and on 



the type of paired bases 0, ^3. ^| • Here we will con- 
sider only models with the later kind of interaction en- 
ergy. In recent years several authors examined this kind 
of models with regard to the thermodynamic behaviour, 
i.e., searching for phase transitions and describing the 
type ofphase involved 0, O, 0, 0; Liu and Bund- 
schuh [3, llJI recently discussed, whether native RNA is 
already in the regime of the thermodynamic limit or finite 
size effect have to be taken into account. In this paper 
we numerically investigate a hybrid model of two well 
known energy models [isl Il6|. i.e., a pair energy models 
where only base pairs are considered regardless of their 
neighbourhood, and a stacking energy model, where only 
consecutive paired bases, i.e., forming a stack, gives an 
energy contribution. It has been claimed that the stack- 
ing energy is more relevant than just the pair energy in 
real RNA . Our model contains terms for both types 
of interactions and allows to move continuously from one 
model to the other. We are interested, whether the two 
limiting models are qualitatively different, in particular, 
whether a phase transition occurs, when moving from one 
model to the other. 

The paper is organized as follows. In section sectionlTll 
we define our model, i.e., we formally define secondary 
structures and introduce our energy model. In section 
IIIBI we explain how to calculate the partition function 
with a dynamic programming algorithm. In section ITTll 
we introduce the observables which we investigate in the 
following section Hvl While in section HV Bl and section 
II V CI we do finite temperature calculations, in section 
lIVDI we use the e-coupling method at zero temperature. 



II. THE MODEL 
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Because RNA molecules are linear chains of bases, 
they can be described as a (quenched) sequence TZ — 
{f'i)i=i,...x of bases e {A, C, G,U}, where L is the 
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length of the sequence. Within this single stranded 
molecule some bases can pair and build a secondary 
structure. Typically Watson-Crick base pairs, i.e., A-U 
and C-G have the strongest affinity to each other, they 
are also called complementary base pairs. Each base can 
be paired at most once. For a given sequence TZ of bases 
the secondary structure can be described by a set S of 
pairs (with the convention 1 < i < j < L), mean- 

ing that bases and rj are paired. For convenience 
of notation we further define a Matrix (S'i j )ij=i....,L 
with Sij = 1 if G S, and Sij = otherwise. 

Two restriction are used: (i) Here we exclude so called 
pseudo-knots, that means, for any {i, j), (i' , j') S S, ei- 
ther i < j < i' < j' or i < i' < j' < j must hold, i.e., we 
follow the notion of pseudo knots being more an element 
of the tertiary structure [l7l |. 

(ii) Between two paired bases a minimum distance is 
required: |j — i| > s is required, granting some flexibility 
of the molecule (here s = 2). 



A. Energy models 

Every secondary structure <S is assigned a certain en- 
ergy E{S); note that this energy in general depends on 
the TZ as well, so it is more precisely to write E{S, TZ), but 
we assume that the structure also includes the informa- 
tion about the sequence. With such an energy model it 
is possible to calculate the canonical partition hmction Z 
of a given sequence TZ by summing over all possible struc- 
tures Z = X]<s but it is computationally more 
efficient to compute it by using the partition functions 
of the subsequences, i.e., by a dynamic programming ap- 
proach. 

Motivated by the observation that the secondary struc- 
ture is due to building of numerous base pairs where ev- 
ery pair of bases is bound by hydrogen bonds, one assigns 
each pair a certain energy e{ri,rj) depending only 
on the kind of bases. The total energy is the sum over 
all pairs 



(1) 



e.g., by choosing e{r,r') — +oo for non-complementary 
bases r, r' pairings of this kind are suppressed. 

Another possible model is to assign an energy Es to a 
pair G 5 iff also (z -I- l,j — 1) G S. This stacking 

energy can be motivated by the fact that a single pairing 
gives some gain in the binding energy, but also reduces 
the entropy of the molecule, because through this ad- 
ditional binding it looses some flexibility. Formally the 
total energy of a structure can be written as 



Es{S) 



E 



Eg Si 



(2) 



assuming that for every pair {i, j) G S the bases and 
rj are complementary bases. The total number t of con- 
secutive base pairs is called the stacking size. Single base 



pairs are not considered as stacks, therefore t > 2 for any 
stack. 

Both types of energy models are discussed in the liter- 
ature [iailflj but, to our knowledge, so far no one has 
discussed a hybrid model. We examine the sum of both 
models at once. 



E{S) -.^ Ep{S) + EsiS) 



(3) 



where the parameters Eg and e{r,r') can be adjusted 
freely, including both models discussed above. Here we 
use 



e(r, r') 



Ep if r and r' are compl. bases 
+00 otherwise 



(4) 



with a pair energy Ep < independent of the kind of 
bases. 

Due to the simple structure of the energy model, e.g., 
the energies depend not on the position of the bases 
within the sequence or whether the paired bases include 
some structure elements like hairpins, the ground state 
is highly degenerated _9, J^J . 



B. Calculating the partition function 

Due to the fact that pseudo knots are excluded from 
our model (see section IH A|l , the calculation of the par- 
tition function can be done recursively. The algorithm 
is similar to that of Nussinov [TsL ITgj . The algorithms 
calculates the elements of two upper-triangular matrices 
iZi,j)i<i<j<L and (^ij)i<i<j<L, where Zij is the parti- 
tion function of the subsequence from base rj to rj under 
the boundary condition that bases ri_i and rj+i arc not 
paired, and Z^.j the partition function under the bound- 
ary condition that bases r^-i and rj+i are complemen- 
tary; Zij is only used as an auxiliary matrix. Then Zij 
can be computed from the partition functions Z and Z 
of smaller subsequences in the following way (remember 
that s denotes the minimum distance between two bases 
of a pair) 

k—i 

Z - Z ■ 1 +e"('=(''-''^)+^=)/'=^^Z-_ui ■ 1 

^t,J ^Z,J — 1 ^ ^ ^t-\-l,j — l 

k=i+l 

Zij = Zij = 1 for i > j 
Zij =0 for < j - i < s - 2 

which is schematically explained in Fig. ^ Because both 
matrices depend on each other, they must be calculated 
simultaneously, starting along the diagonal and continu- 
ing along the off-diagonals. The calculation of the par- 
tition function can be done in 0{L^) steps, where L is 



3 



E 



(i. k - 1) 



k + l,j-l) 



{i.i ~ 1) 



(i,fc - 1) 



k + l.J-l) 



--E„ + E,, 



FIG. 1: Schematical explanation of eq. Q and eq. JSJ, e.g., 
white boxes represent Z, gray boxes Z. 



the length of the sequence. The partition function Z of 
the entire sequence is Zij^, but also the other matrix el- 
ements are useful for generating ensembles of structures 
according to the Boltzmann distribution (see section llTT|l . 

A similar algorithm can be used to calculate the ground 
state energy: 



Ni^j = min(7Vi.j_i, min {Ni,k-i + e{rk,rj)) + Nk+i.j-i) 

k—i 

j-s-i . (6) 

min {Ni^k-i + e{rk,rj) + Nk+ij^i)) 

k=i+l 

Ni^j = Ni,.j =0 for i > j 

J = +00 for < j - i < s - 2 

In comparison to eq. |SJ| additions are replace by min- 
operations, multiplications by additions and the expo- 
nentials of the energies by the energies them self. 



III. OBSERVABLES 

After calculating the partition function Z for a given 
random sequence, we want to measure some quantities 
to compare the members of the ensemble. In principle 
most quantities could be calculated by a similar dynamic 
programming algorithm introduced above, but in gen- 
eral the running time behaviour would be of higher order 
(than three) in the sequence length. For this reason we 
use a different method, where an ensemble of structures 
is generated due to its Boltzmann weight The pro- 
cedure to build a sequence is essentially a backtracking 
algorithm: Starting with the entire sequence a partner 
for an outermost base, e.g., base L, is chosen with the 
appropriate probability, e.g., base fc, after this the proce- 
dure is applied to the subsequences 1 to A: — 1 and fc -|- I 
to L. If base L is chosen not to be paired at all, one 
uses the sequence shortened by base L. Considering a 
subsequence from base k to Z, the probability pij:^^; that 



bases i and j {k < i < j < I) are paired is given by 

( Z-^ Zi ■ ,e'i<^''^-'^)+^°)/k^T ^ 

k.l — + 

bases i ~ 1 and j + 1 paired 



Pi,j;k,l 



bases i — I and j + 1 unpaired 



(7) 



For each member of this ensemble £ the quantity 
of interest X is calculated and the average {X) — 
j^^^s X{<S) is used as an approximation to the expec- 
tation value of the Gibbs-Boltzmann-ensemble; for large 
enough ensembles £ this average approaches the true ex- 
pectation value. 

A simple observable is the energy E and its fluctuations 
(Ai5)^, the later is directly connected with the specific 
heat cv = {AEf/LkBT^- 

Of particular importance is the overlap q between two 
structures S and S' 



q{S, S') := 



- y 

l<i<j<L 



(8) 



that is the number of bases paired to the same base in 
both structures normalized such that q lies always be- 
tween zero and one. This is a measure of how similar 
two structures are. Note however, that with this defini- 
tion the overlap of one structure with itself is q{S, 5) ^ 1 
unless all bases are paired, where it is equal to one. A 
definition of q where also bases unpaired in both struc- 
tures are counted is used in resulting in an over- 
lap definition that is normalized, however this similarity 
measurement has the drawback that the less bases are 
paired the more two structures get similar. We further 
remark that for any two structures S, S' the Cauchy- 
Schwarz inequality (g(5,5'))^ < q{S,S) q{S' ,S') holds. 
With this quantity q two ensembles £ and £' can be com- 
pared, e.g., looking at the distribution of q{S, S') of all 
Se£,S' e £'. 

The ensemble averages in general depend on the 
chosen sequence, therefore a further averaging over sev- 
eral (random) sequences is required. This sequence aver- 
age is denoted by [•]. We again approximate this average 
by summing over a finite set of sequences. 

Because of this two stage averaging, it is probably 
preferable instead of looking at [(g)] directly to use func- 
tions of the first and higher moments of 9, e.g., the Binder 
cumulant pol I21I , which is defined by 



(9) 



where (g") is either the average over all pairs of one en- 
semble 




\£\(\£ 



(10) 



s, s'e£ 
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or the average over all pairs of two ensembles 



Basic observables 



\e\\£ 



se£,S'e£' 



The later choice is appropriate when one is looking for a 
change in the behaviour of the models when one varies 
the parameters in comparison to a reference model, while 
the former is the better choice, if an external parameter, 
e.g., the temperature, is varied. 

The Binder cumulant B vanishes at high temperatures, 
while for low temperatures it approaches a finite value in 
the thermodynamic limit. 

A similar quantity has been used in 0: 



(11) In section lll Al where we introduced the energy model, 

we opposed the pair energy to the stacking energy model. 
In Fig. 13 we show the average size t of a stacking as 
a function of the energy parameter Ep at temperature 
T = 0. We keep Ep + Eg = —1 constant, to fix the 
overall energy scale. For all fixed energy parameters the 



A := 



[xV\ - [xn] 
[Xnf 



where 



xn -.^ L [{q^) - {qf) 



(12) 3 



(13) 



is the variance of the q distribution. This parameter A 
measures how the probability distribution of q varies from 
sequence to sequence. A value close to zero indicates a 
self- averaging behaviour. 



IV. NUMERICAL RESULTS 




Average stacking size as a function of energy param- 
jp for different system sizes. At Ep = —1, i.e., Es = 0, 



FIG. 2: 

eter Ep 

and Ep = the curves are discontinuous. For explanation of 
symbols see caption of Fig. |7| Missing error bars are of the 
size of the symbols or smaller, and omitted for legibility. 



In order to find some possible differences in the be- 
haviour of the energy model eq. Q for different param- 
eters Ep, Eg, as introduced in eq. Q and eq. (0), re- 
spectively, we numerically calculated the quantities men- 
tioned in section IIIII above. In all our examples we av- 
eraged over randomly generated sequences TZ = (r^), 
where the probability for a specific base at position 
i is J for all base types independent of the other bases 
Tj^i. The size of the sequence varied from L = 128 up 
to L = 1024, for the disorder average 2000 up to 6000 
random sequences were used. Pairing of bases are only 
allowed for complementary bases and the minimum dis- 
tance between to bases was chosen as s = 2. 

In section lTV Al and If V lj| we vary continuously the en- 
ergy parameters between the two extreme cases (Ep — 
0,Es = -1) and {Ep = -l,^;^ = 0). In section 117^ we 
shortly discuss the averages of the stacking size t and of 
the q for different energy parameters. In section HV Iji we 
examine the cross overlap qg gici, see eq. (|ll|l . between 
a reference ensemble S'^'^^ generated for fixed energy pa- 
rameters, and ensembles generated for different energy 
parameters. In the following section HV CI we look at the 
temperature variation of various quantities without using 
any reference ensemble to estimate some critical param- 
eters. In the last section HV Dl we apply the e-coupling 
method at T = to determine the critical exponent 9 
describing the behavior of low-lying excitations. 



average stacking size t increases with increasing system 
size, which is expected as with increasing length the prob- 
ability for longer stacks increases. Also as expected is the 
overall increase in the average stacking size with the en- 
ergy parameter Eg, because the stacking energy prefers 
to build stacks. The large increase of t while changing 
Es from zero to a nonzero value can be explained as fol- 
lowing: The ground state for Ep = —1.0, Eg = 0.0 is 
highly degenerated, while changing to a nonzero Eg only 
those states stay ground states which have a high stack- 
ing contribution to the energy. This selection increases 
the average of t. A similar argument applies at the other 
end, where Ep changes from a nonzero value to zero. 

Similarly the overlap q jumps to a larger value when 
changing from Ep = to a nonzero value (Fig. In 
addition at positions, where Ep/Eg are fractions with 
small numerator and denominator, e.g., y or i, for this 
energy parameters the ground states in configuration 
space are more broadly distributed than for slightly dif- 
ferent parameters. This can be seen in the right inset 
of Fig. 13 where g-Distribution in the symmetric case 
{Eg = Ep = —0.5) is broader than in the slightly asym- 
metric case {Eg = -0.49, Ep = -0.51). The width of this 
minima in the main plot, as well as that in Fig. ^ and 
Fig- El is below AE = 0.001, as one can estimate from 
the left inset of Fig. El 

For both the average stacking size and the average 
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FIG. 3: Average overlap qs as a function of the energy pa- 
rameter Ep for different system sizes. The local minima at 
Ep = —0.75, —0.5, —0.25 are due to the commensurability of 
Ep and Ea and indicate a broad distribution of the ground 
states in configuration space. For explanation of symbols see 
caption of Fig. |7| The left inset is an enlargement to show 
the discontinuity. The right inset is the g-distribution for 
Es = Ep = -0.5 (solid lines) and Es = -0.49, Ep = -0.51 
(broken lines) for sequence lengths L — 512 (thin lines) and 
L = 1024 (thick lines). 



overlap, the behavior changes smoothly when moving 
from one model to the other, with the exception of the 
highly degenerate points, where we can observe the jumps 
in the overlap q. Hence, there is no sign of a transition 
in between the two extremal models. To confirm this, we 
next study the Binder cumulant. 



B. Binder cumulant 

Since we introduced in eq. lj3Jl a whole class of energy 
models depending on the pair energy Ep and the stack- 
ing energy Eg , we examined the behaviour of the Binder 
cumulant depending on this two energy parameters and 
the sequence length L. Second order phase transitions 
are characterized by crossing of the Binder cumulant for 
different system sizes at the transition point. 

In Fig. 2] the Binder cumulant cq. Q is shown at T = 
using the "self-overlap" eq. (|10|l . Again, the energy 
parameter Ep and E^ are varied such that always Ep + 
Es = —1 holds. The value of B increases with increasing 
system size for Ep, i.e., the curves do not cross each other 
and therefore give no evidence of a phase transition. The 
local minima are - as the minima in Fig.|31- due to the 
commensurability of the energy parameter Ep and E^. 

Further we used a reference ensemble generated for 
energy parameters Ep — Es — —0.5 and used the "cross 
overlap" of eq. Hll|) . First, the values of B at Ep = Eg — 
—0.5 coincide with the values shown in Fig. 0] Similar 
to the observation above, B roughly increases with the 



FIG. 4: The Binder Cumulant B of eq. ® with q ^ qe 
(see eq. llOl 'l at temperature T = 0. The curves for different 
system sizes do not cross, and therefore B gives no hint for 
a phase transition. For explanation of symbols see caption of 
Fig. El 
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FIG. 5: The Binder Cumulant B of eq. @ with q = qe,e' (see 
eq. moi 'l at temperature T = 0. For explanation of symbols 
see caption of Fig. |7| 



system size, although the separation of the curves is not 
as clear as in Fig. ^ especially in the range from Ep = 
—0.4 to Ep = 0.0, where the curves coincide within the 
error-bars. To summarize, the dependence of the Binder 
cumulant on the energy parameters does not indicate a 
phase transition. 



C. Temperature dependence of the energy models 

Another possible method to discriminate between the 
different energy parameters is to examine the tempera- 
ture dependence of some quantities, especially the be- 
haviour at a critical temperature. In Ref.i it was shown 
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for a similar model that below a critical temperature, al- 
most all sequences fold to a compact structure, but for 
most sequences not into a single structure. They point 
out that this kind of low temperature behaviour is well 
known from spin glass and other disordered systems. In 




T 



FIG. 6: Specific heat cv as a function of temperature for 
parameters Es = —1.0 and Ep = 0.0. The curves for different 
system sizes crosses at T « 0.25 and T ~ 0.11. The inset is 
an enlargement of the main plot. For explanation of symbols 
see caption of Fig. |7| 

Fig. we plotted the specific heat for different system 
sizes as a function of temperature for Eg ~ —1.0 and 
Ep = 0.0. All curves cross each other close to T — 0.11 
and T — 0.25, which is an evidence for a phase transition 
at this temperature region. The data for other energy 
parameters {Es,Ep) look similar, but the curves cross at 
different temperatures. 

To determine the critical behaviour quantitatively we 
investigated the width xn of the overlap distribution. As 
can be seen from Fig. \7\ all curves have a maximum and 
the position of this maximum is decreasing with increas- 
ing sequence length. We assume that the maximum po- 
sition has the form Tc{L) = Tc + a L^^^'^ and fit the data 
to this form to get the critical parameters (see Fig. |HJ). 
The results for three different pairs of energy parameters 
are shown in Tab. ^ The critical exponent v for the pa- 



Energy Model 


\lv 




e 


Es=0 
Es = -0.5 
Es = -1 


Ep = -1 
Ep = -0.5 
Ep = 


0.93(15) 
0.70(36) 
0.36(17) 


0.086(3) 
0.109(7) 
0.125(21) 


0.229(38) 
0.237(50) 
0.194(67) 



TABLET: Critical parameter of a XTC-maximum fit. Compar- 
ing the two limiting cases {Es,Ep) = (—1,0) and {Es,Ep) — 
(0,-1) the parameters v and Tc are different and indicate 
a quantitative change. The last column belongs to the e- 
coupling method in section HV Dl 



rameter pair {Es ~ 0,Ep = —1) is clearly different from 
that of the parameter pair [E^ = —l,Ep = 0), showing 
that the quantitative behaviour changes when varying 




Q\ , \ , \ , 1 , \ , 1 , 1 

0.08 0.1 0.12 0.14 0.16 0.18 0.2 

T 

FIG. 7: Xn ^ ^ function of temperature for energy parameter 
Es — —1.0, Ep = 0.0. For system size L the following symbols 
are used: o 128, A 192, □ 256, V 384, 512, <s 768, ■ 1024. 
Calculated data points are indicated by symbols. Lines are 
drawn to guide the eye. Missing error bars are of the size of 
the symbols or smaller, and omitted for legibility. 
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FIG. 8: The position of the maxima of X'R- for differ- 
ent energy parameter. The curves are fitted to the form 
Tc{L) = Tc + aL~^^" , and gives the critical parameters of 
Tab. H] Missing error bars are of the size of the symbols or 
smaller, and omitted for legibility. 



from one limiting case to the other. On the other hand 
for the intermediate parameter set the critical exponent 
is consistent with both within the error range. 

Finally, we also found (not shown) that the behaviour 
of A of eq. p2(l is in agreement with this observation and 
with the results for a two-letter RNA model. For all 
three cases studied here, the width of the g-distribution 
varies only slightly from realization to realization at high 
temperatures, while for low temperatures the self aver- 
aging behaviour disappears. 
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D. e-coupling 

Previously the e- coupling method has been used |l2ll23| 
to investigate low-energy excitations of RNA secondary 
structures. The basic idea is to add another term to the 
energy model in eq. which depends on the ground 
state of the original problem: Assume Sq is the unique 
ground state of E{S), then a new energy function is de- 
fined as following 

E,iS)^EiS)+eqiS,So) (14) 

with e > 0. The additional term penalizes structures 
similar to So, where eq{S,So) is largest for S = So- In 
general the ground state S^ of the new energy model 
will be different from Sq. The difference AE{e) := 
E{S^) — E{So) is an increasing function of e and A£^(e) < 
£ holds. The later implies that S^ is for small enough e 
a low lying excitation of the original energy model, and 
has the smallest overlap with Sq- 

The average distance d{e,L) := 1 — q{Se,So) between 
the new and the old ground state scales as d{e,L) ex 
e , £ constant, while the average energy difference 
scales as AE{d,L) cx L^, d constant, with the critical 
exponents 9. For details see Refs. IT^ orl2^ 

The assumption of a unique ground state used above 
does not hold for our model used so far: in general the 
ground state is highly degenerated because only two en- 
ergy parameters (Es and Ep) are used, many structures 
will have the same energy. The degeneracy of the ground 
state renders the e-coupling method as described above 
almost useless. Since in natural RNA the contributions 
to the energy are more complex different structure will 
never be degenerated. This justifies to change the energy 
model slightly by adding a random energy r]ij to the pair 
energies introduced in eq. Q): e{ri,rj) e{ri, Vj) -f rji^j. 
There are several possibilities to choose the distributions 
of the r]s (see 12]), here we use identical independently 
distributed Gaussian random number with zero mean 
(t/) = and variance (rf) = r^^jL with 770 = 0.1 (the 
QD model in Ref. 12). 

The raw result for different values of e G [0.01, 100] is 
shown in Fig. |^ A scaling plot of the data for e < 1 
according the scaling form A.EL~^ = f{d) is shown in 
the inset of Fig. El The scaling parameters leading 
to the best data collapse for different energy parameters 
are shown in the right most column of Tab. ^ They are 
equal within the error margins and thus does not give 
us a further hint of a quantitative different behaviour for 
different energy parameters. This difficulties in doing a 
good scaling of the data in the QD model were already 
pointed out 0. However, for a different model using 
a scaling function with finite-size corrections a critical 
exponent 9 = 0.23 ± 0.05 was obtained |2^, which is 
close to our findings. 
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FIG. 9: e-coupling results for Es = -0.5, Ep = -0.5. The 
energy difference AE(S) between the original ground state 
structure iSo and the ground state structure of the dis- 
turbed model is plotted as a function of the distance between 
this structures. The inset is a scaled plot of the data with 
£ < 1 of the main plot. [6 = 0.24, see Tab.HJ For explanation 
of symbols see caption of Fig. 

V. SUMMARY 

We have introduced a RNA model which continuously 
interpolates between two well known models. We sought 
for the answer to the question, whether there is any phase 
transition of the thermodynamic behaviour. We used 
both zero as well as finite temperature data. 

Zero temperature results give no evidence for a phase 
transition apart from trivial transitions, e.g., discontinu- 
ity of t at points with Eg ^ or Ep — Q. The curves 
of the Binder cumulant do not cross at a certain point, 
which would be an indication of a phase transition. Sim- 
ilar, the critical exponent 9 derived from the £-coupling 
method seems to be independent of the energy param- 
eters Es and Ep, and therefore gives no evidence for a 
quantitative difference in the thermodynamic properties. 
But as stated in ^3 , the determination of the critical pa- 
rameter is rather difficult in the quasi-degenerated case 
studied here. 

The finite temperature results show - in contrast to 
the zero temperature data - a quantitative dependence 
on the energy parameters. The critical exponent v for the 
correlation length seems to depend on the energy model 
and may vary continuously while going from one limiting 
model to the other. 
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